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ABSTRACT 



We discuss images of the central ~10 kpc (in projection) of the galaxy merger 
NGC 6240 at H and K' bands, taken with the NIRC2 narrow camera on Keck 
II using natural guide star adaptive optics. We detect 28 star clusters in the 
NIRC2 images, of which only 7 can be seen in the similar-spatial-resolution, 
archival WFPC2 Planetary Camera data at either B or I bands. Combining the 
NIRC2 narrow camera pointings with wider NICMOS NIC2 images taken with 
the F110W, F160W, and F222M filters, we identify a total of 32 clusters that are 
detected in at least one of these 5 infrared (A c > 1 /mi) bandpa sses. By comparing 



to in stantaneous burst, stellar population synthesis models (IBruzual fc Chariot 



20031 ). we estimate that most of the clusters are consistent with being ~15 Myr 
old and have photometric masses ranging from 7 x 1O 5 M to 4 x 1O 7 M . The total 
contribution to the star formation rate (SFR) from these clusters is approximately 
1OM yr _1 , or ~10% of the total SFR in the nuclear region. We use these newly 
discovered clusters to estimate the extinction toward NGC 6240's double nuclei, 
and find values of Ay as high as 14 magnitudes along some sightlines, with an 
average extinction of Ay ~ 7 mag toward sightlines within ~ 3"of the double 
nuclei. 
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INTRODUCTION 



Ultraluminous Infrared Galaxies (ULIRGs) output more energy in the infrared than at 
all other wavelengths combined, and have Ljr > 1O 12 L . In the local universe, these relatively 
rare objects are frequ ently associated with peculiar galaxies, thought to be merging systems 
(ISanders et al. Ill988l ). Their high infrared emission comes from warm dust, heated by a 
combination of starbursts and active galactic nuclei (AGN) which are thought to be fueled 
by huge quantities of molecula r gas that fall inward tow ard the nucleus during a major merger 
of two gas-rich spiral galaxies (ISanders fc Mirabellll996l). While at the pre sent epoch ULIRGs 
are rare in number, numerical simulations by iMihos fc Hernquistl ( 119941 ) suggest that they 
are not merely anomalies resulting from "special" merger conditions, but rather their rarity 
is due to the fact that they represent a short-lived yet common phase of galactic evolution, 
namely when spiral galaxies collide to form massive elliptical galaxies. Detections by the 
SCUBA telescope o f hundreds of high -redshift sub-mm galaxies (SMGs) lend further support 



to this th eory (e.g. Ifjorvs et al. 



20031 ). The properties of SMGs are similar to those of local 



ULIRGs (jConselice et al. 1120031 ). and the number density of SMGs is similar to the number 
density of massive elliptical galaxies at the present epoch. Thus, assuming that structure 
in our Universe forms by hierarchical assembly, ULIRGs offer us the unique opportunity to 
glimpse a brief but important stage in a galaxy's life, when massive spirals transform into 
ellipticals. 

Though NGC 6240 (z = 0.024, d = 98 Mpc, 1" = 477 pc; H = 71 km s" 1 Mpc" 1 , 
Qm = 0.27, Qa = 0.73) technically lies at the high-luminosity edge of the LIRG class with 
Lir ~ 6 x 10 11 Lq, in many ways it is a prototypical ULIRG, and is commonly referred to as 



such. NGC 6240 is a late-stage merger of two massive disk galaxies ( ITacconi et al. Ill999l ); it 
has extended t idal tails and a close nuclear pair separated 



2001 



of 715 pc (e.g. Max et al. M20051 . 120061 : iBeswick et al. 
by red supergiants in ~20 Myr-old starbursts (ITecza et al. 



by ~1.5", or a projected distance 
). Its K' band light is dominated 



2000) , and at the double nucleus 



two AGN have been discove red in hard X-rays (IKomossa et al. 1120031 ) as well as at 5 GHz 
(IGallimore fc Beswickl|2004l). As in other mergi ng systems, most notably NGC 4038/9 (the 
Antennae galaxy) (IWhitmore fc Schweizerlll995l ). star clusters have previously been observed 
in NG C 6240's main body a nd tidal tails with the WFPC2 F 450W, F547W, and F814W 
filters (IPasquali et al. II2003I ). But recently iMax et al. I (120051 ) detected previously unseen 
circumnuclear star clusters at J, H and K' bands using natural guide star (NGS) adaptive 
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optics (AO) with the NIRC2 wide-field camera on Keck II. These clusters, and other similar 
clusters in ULIRGs and galaxy mergers, may hold clues regarding a variety of open questions, 
from the distribution of dust and star formation in ULIRGs, to the formation and evolution 
of globular clusters (GCs). 

Historically, measuring star formation rates (SFRs) in galaxy mergers has been difficult, 
and until recently the contribution of SMGs to the total luminous energy density in the 
high-redshift universe has been significantly u nderestimated. A variety of heterogeneous 
techniques have been used to measure SFRs fcennicuttl[l~998L but in starburst galaxies 
these are complicated by the presence of dust (jCalzettil l200ll ) . including the usual methods 
of measuring SFRs through Ha line fluxes or UV continuum luminosities. In these cases 
astronomers have turned to far-infrared (FIR) wavelengths to measure the dust-reprocessed 
light. But a well-known problem arises when using FIR luminosities to estimate SFRs in 
ULIRGs: in addition to young OB stars, AGN also help heat the dust. While studies 
have shown that the lu minosities of most ULIRGs are dominat ed by star formation even 
when AGN are present ( iGenzel et al. 1 1 19981 ; iKewley et al. 1 120011 ) , accurately measuring the 
fractional AGN contribution to the bolometric luminosity (Lboi) has proved difficult. For 
example, lArmus et al.l (120061 ) reported on their detection of the [Ne v] 14.3 /im emission 



line feature, with a flux of 5 x 10~ 14 erg cm -2 s~ l , which provides direct evidence of the 
buried AGN in NGC 6240. They note that the [Ne v]/[Ne 11] ratio is often incorrectly 
treated as a reddening independent quantity, even though [Ne v] comes directly from gas 
heated b y the AGN, while [Ne 11] can be dominated by a starburst. Without correcting 
for dust, lArmus et al. state that the large equivalent width of the 6.2 /im polycyclic 

aromatic hydrocarbon (PAH) emission feature, and the small [Ne v]/[Ne 11] and [Ne v]/IR 
flux ratios, are consistent with an apparent AGN contribution of only 3 — 5% of Lboi- But 
after c orrecting the measur ed [Ne v] flux for extinction using values implied by their SED 
fitting, lArmus et al.l (120061 ) estimate that the intrinsic fractional AGN contribution to the 
bolometric luminosity of NGC 6240 is much higher, between 20 and 24%. 

By detecting individual, massive, star clusters in merging galaxies we can trace the most 
violent star formation activity, while circumventing some of these hurdles which make accu- 
rately measuring SFRs a difficult task. For example, the circumnuclear locations of clusters 
makes them easily distinguishable from AGN. In addition, while the exact size distribution 
of dust grains in both star clusters and diffuse stellar populations is impossible to know, the 
complicated dust mixing models required for analyzing diffuse stellar fields reduce to simple 
dust screens toward clusters, if we assume that all intracluster dust has been evacuated. 
Even when some intracluster dust is present, the dust screen approximation is still reason- 
able assuming that the amount of dust internal to the cluster is far less than the intervening 
dust along the line of sight through the entire merging galaxy. Furthermore, individual clus- 
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ters are well approximated by single stellar populations, making stellar population analysis 
of clusters more straightforward than for diffuse stellar fields. Finally, near-infrared (NIR) 
cluster observations provide a more localized measure of the SFR compared to estimates 
based on longer wavelength observations. But unfortunately, while star clu sters seem to 



represent a common mode of star f ormation in spirals and LIRGs alike (e.g. iLarsenl 12002 



Scoville et al. Il2000l ; iFall et al.ll2005l ) , detecting dust-reddened clusters is challenging because 



longer wavelength observations are innately of lower spatial resolution. However, with the 
improvement of AO systems on large, ground based telescopes, the possibility of detecting 
individual infrared star clusters, even at z ~ .02, is now a reality. 

In order to measure localized SFRg^ with fewer complications due to AGN and dust, we 
have obtained ~ .04" resolution images of roughly 10 kpc surrounding the merging galaxies' 
nuclei using the NIRC2 narrow-field camera with NGS AO at Keck. We combine our data 
with HST images taken with the NICMOS NIC2 camera (using the F110W, F160W, and 
F222M filters) and the WFPC2 Planetary Camera (using the F450W and F814W filters). 
Using the combined data set we study a total of 32 clusters with inhomogeneous wavelength 
coverage, partially due to variations in the fields of view, and partially due to non-detections 
in some wavebands. In order to estimate the SFR associated with these clusters we es- 
timate masses, ages, and extinctions toward these 32 clusters by comparing the clusters' 
spectral energy distributio ns (SEDs) to model SEDs from stellar population synthes i s code s 
(IBruzual fc Charlotl 120031 ) . employing simple, screen extinction models by iDraind (120031 ). 
Given the data currently available, we do not attempt to quantify the cluster's properties 
with the precision of previous studies of nearby mergers. This work is meant to be a prelimi- 
nary analysis of NGC 6240's recently discovered cluster population - one of the most distant 
such populations ever detected; future high-resolution spectral observations will provide a 
much clearer understanding of the cluster's properties. In §2] we describe our observations, 
in §3] we describe our data reduction, including the techniques we used to identify clusters 
and perform photometric calibrations, and in §H and we present and discuss our results. 



1 We note that this method of using star clusters to measure localized SFRs, while beneficial in that it 
sidesteps complications due to AGN and dust, by no means traces the total SFR. We preferentially trace 
star formation in the most massive, least dusty environments. Less massive or more obscured clusters likely 
remain undetected in this work. 
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2. OBSERVATIONS 

2.1. Adaptive Optics at Keck 

We observed NGC 6240 with the NGS AO system on the Keck II telescope using the 
narrow-field camera (0.009942" pixel -1 ; 4.74 pc pixel -1 ) of the NIRC2 infrared instrument at 
H (A c = 1.633/im) and K' (A c = 2.124/xm) bands. The H and K' observations were obtained 
on 24 April 2005 and 17 August 2003, respectively. The natural guide star (B = 13.5) is 
USNO-B 1.0 0924-0386013, located 35.8" east-northeast of NGC 6240's dual nuclei. 

Due to the faint nature of the star clusters superposed on the diffuse background light 
from the galaxy, it is problematic to measure the resolution and Strehl ratio (S) at the 
location of the clusters, i.e. off-axis from the guide star. We did not observe PSF star pairs 
and therefore cannot reasonably guess what the Strehl may be at the clusters' locations. If 
there were no tip-tilt errors then the core of the off-axis point spread function (PSF) would be 
diffraction-limited, and at 2.1/im the full-width-half-maximum (FWHM) of this diffraction- 
limited, off-axis PSF core is ~ 0.04" (~ 19 pc at z = .024), making our observations 
slightly oversampled. Throughout this paper we assume that the off-axis PSF and Strehl 
- whatever they may be - are constant across the off-axis FOV. This assumption is well 
justified when a small FOV (~ 10") is far off-axis (>30"), and when the areas of interest in 
the FOV are oriented perpendicular to the direction of the guide star. Most of the clusters 
in our images are located along the north-south orientated double nuclei, and therefore 
are oriented approximately perpendicular to the direction of the guide star. The assumed 
constancy of the resolution and Strehl ratio across the FOV implies that our detection limits 
and photometric measurements should be accurate in the relative sense; in §3.41 (case II) we 
describe the steps we have taken to ensure that our photometric measurements are accurate 
in the absolute sense as well, by sidestepping the need to know the absolute Strehl ratio to 
perform photometry. 

In Table [1] we outline the details of these Keck observations, including the on-axis Strehls 
achieved for the guide star on each night. While the off-axis Strehl is unknown, the Strehl 
measured on the guide star is an indicator of the overall performance of the AO system, 
including the effects of wind turbulence, sampling rate, and weather. 

2.2. Archival HST Observations 

We have analyzed archival WFPC2 images taken as part of the GO proposal 6430 (PI: 
R. van der Marel). The observations were taken with the Planetary Camera (.0455" pixel -1 ; 
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21.7 pc pixel" 1 ) using the F450W and F814W filters on 28 March 1999. The theoretical 
FWHM of the diffraction-limited PSFs for these two WFPC2 images are ~ .038" and ~ .069", 
respectively. At NGC 6240's distance, these angular resolutions correspond to 18.1 and 32.7 
pc, respectively. The F450W and F814W filters are roughly equivalent to the Johnson- 
Cousins B and I filters, therefore throughout this paper we use the shorter Johnson- Cousins 
notation interchangeably with the true WFPC2 filter names. However, when performing all 
calculations and measurements we use the exact filter transmission functions for the WFPC2 
filters. 

We have also made use of archival NICMOS images obtained on 1998 February 02 
with the Hubble Space Telescope using the Near Infrared Camera and Multi-Object Spec- 
trometer i _(NICMOS)_as i_p_art of the Guaranteed Time Observing (GTO) program 7219 (PI 



Scoville; IScoville et al. Il2000l ). The imaging program employed three spectral filters (F110W 
A c = 1.100 /on, FWHM = 0.592 /im; F160W: A c = 1.594 /im, FWHM = 0.403 /im; F222M 
A c = 2.116 /im, FWHM = 0.143 /im) in NICMOS camera 2 (pixel scale: x sca i e = -0762" 
pixel" 1 , y S caic = .0755" pixel" 1 ; x sca i e = 36.3 pc pixel" 1 , y sc& i e = 36.0 pc pixel" 1 ). Four 
exposures in each filte r, using STEP8/NSAMP = 10, 11, and 12 multiaccum sampling 



(IMacKentv et al. 1 119971 ) for F110W, F160W and F222M, respectively, were taken with cor- 



responding integration times of 40s, 48s, and 56s in eac h exposure. NIC MOS camera 2 



critically samples its point spread function at 1.6 /im (see iMax et al. 1 120051 . §2.2). All im- 
ages were acquired under a four-point "dither" pattern with pointing offsets of 1.9215", thus 
enabling critical sampling at 1.1 /im (after image combination), as well as improved sam- 
pling and rejection of data from defective pixels at all wavelengths. Because the F160W 
and F222M filters provide wavelength coverage roughly similar to the H and K' Mauna Kea 
filters (Fig. [1]), these two NICMOS images are used primarily as a tool to help perform pho- 
tometric calibrations on the NIRC2, AO images. (The details of this photometric calibration 
technique are described in §3.41 ) However where the fields of view of the NICMOS F160W 
and F222M images extend beyond the NIRC2 fields of view, the NICMOS images provide 
additional spectral information and thus we analyze these images. Note that of all the data 
discussed in this paper the NICMOS images have the lowest spatial resolutions - about 4 
times lower than the resolution provided by the Keck AO system at similar wavebands. The 
theoretical FWHM of the diffraction-limited PSFs for the three NICMOS images are 0.093", 
0.13", and 0.18", corresponding to linear scales of 44.2, 64.3, and 85.0 pc. 

Details regarding both the WFPC2 and NICMOS HST observations are shown in Ta- 
ble m 
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3. DATA PROCESSING 

The complete set of data is composed of 7 images from 3 different instruments at the 
following bandpasses, in order of increasing central wavelength: B, I, F110W, F160W, H, 
K', and F222M. The inhomogeneous nature of these images, including the manner in which 
they were taken, the morphology that they display, and their resolutions which span more 
than a factor of five, requires that we employ non-uniform processing techniques. We shall 
discuss the pitfalls inherent to treating data in a non-uniform fashion. 



3.1. Data Reduction 

The near-infrared images taken with NIRC2 were reduced in the standard way, by 
masking bad pixels, flat fielding each frame using twilight images, and dark subtracting and 
sky subtracting frames of the appropriate exposure times. Before combining NIRC2 narrow 
camera frames, we removed the geometric distortions using the IRAF DRIZZLE package 
with the distortion coefficients published in the NIRC2 manual, and the kernel parameter 
set to "lanczos3." 

The NICMOS data analyzed in this paper were obtained by recalibrating and repro- 
cessing the archival data. Instrumentally calibrated count-rate images were created from the 
raw archival multiaccum frames after applying linearity corrections, dark-frame subtraction, 
DC bias compensation, flat-field correction, and cosmic-ray rejection. Darks, made from 
contemporaneous on-orbit calibration observations, calibration reference files for linearity 
corrections, obtained from STScI, and flats, obtained from the NICMOS Instrument Defini- 
tion Team (IDT), were applied in an IDL-based process closely following the IRAF/STSDAS 
calnica task. The image data were photometrically calibrated based upon the absolute in- 
strumental calibration established by the IDT (F110W: 2.031 x 1(T 6 Jy count -1 s; F160W: 
2.190 x 10 -6 Jy count" 1 s; F222M: 5.487 x 10 -6 Jy count" 1 s). Known "bad" or defective 
pixels were replaced by two-dimensional Gaussian weighted interpolation of good neighbor 
pixels with wavelength-dependent weighting radii of 3, 5, and 7 pixels for the F110W, F160W, 
and F222M frames respectively. The fully calibrated images were geometrically rectified to 
correct for the ~ 0.9% linear geometrical distortion in NICMOS camera 2, and for each filter, 
those four images were median combined into a single image on a 2-times finer resampled 
grid (~ 38 mas pix" 1 ). All image resampling (including re-mapping for registration, and 
rectification) was done with sinc-function apodized Gaussian weighted interpo lation using 



the IDP3 package developed by the NICMOS IDT ({Schneider fc Stobid 120021 ). conserving 
the target flux density while resampling. The three resulting images (one for each filter) are 
those that are discussed further in this paper. 
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We obtained the archival WFPC2 data from the WFPC2 associations at MAST (Multi- 
mission Archive at STScI), and these data are corrected for bias, dark current, flat fielding, 
and cosmic rays. We performed no additional processing on these WFPC2 associations. 



3.2. Cluster Identification 

We attempted to identify unresolved star clusters in the 7 images in a systematic way 
by using the FIND routine adapted for IDL from DAOPHOT. However, this method and 
variations on this method, including using the FIND algorithm on unsharp-masked images, 
were all unsuccessful, especially in regions where the diffuse background light from the galaxy 
is highly variable or bright. Therefore in this paper we have visually identified clusters, and 
we define a cluster to be any source that appears unresolved in an unsharp masked image. We 
neglect clus ters that are visible on ly at B and/or I bands because these have been previously 



analyzed by lPasquali et al. I (120031 ). and our analysis would provide no new spectral coverage. 



Figure [2] shows the cluster positions superimposed on images at B, I, and K' bands. 

In total we identified 32 clusters that are visible in at least one of the infrared passbands 
(F110W, F160W, H, K', and F222M). Of these, 27 lie within the intersection of the 7 fields 
of view, lending themselves to the most complete spectral analysis. Of these 27 clusters, only 
7 have low enough extinction to be detected at I band, and this number decreases to 5 at B 
band. Figure [3] shows the locations of all 32 clusters and the relative sizes and orientations 
of the 7 fields of view. In Figure H] we label each cluster with a number from 1 through 32 
with 1 being the northernmost cluster. N ote that clus t er #1 5 is coincident with the site of 



the northern black hole, as determined by iMax et al. I (120061 ). For reference, in the publicly 



available WFPC2 I band image, the brightest pixel in the northern nucleus is 0.27" east 
and 0.099" north of the northern black hole (i.e. cluster #15). In Figures |3] and H] we mark 
the position of the northern black hole with an orange cross. The positions of each cluster 
relative to cluster #15 can be found in Table [31 

The complicated morphology of NGC 6240 makes it impossible to specify a single cluster 
detection limit across the entire FOV. The morphology of the H and K' images is relatively 
smooth, with background luminosity increasing sharply toward the double nuclei. Therefore 
the cluster detection limits at H and K' worsen monotonically towards the double nuclei. 
While it is generally true that the detection limits worsen toward the double nuclei in the 
B and I images as well, the situation here is more complicated because these images are far 
more patchy due to dust obscuration. We estimate our detection limits by adding simulated 
clusters to various locations in each image, and decreasing the fluxes of the simulated clusters 
until the clusters become undetectable by eye. (See §3.41 Cases I and II for a detailed 
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description of how we create the simulated clusters.) Using this method, we estimate that 
clusters brighter than K\ B = 22 mag, H AB = 23 mag, Jab = 24 mag, and B AB = 26 mag 
should be detectable throughout much of the images. 

3.3. Image Alignment 

In order to align all 7 images to the same spatial grid we first visually determined which 
star clusters are common to pairs of images. (No stars are visible in the small NIRC2 FOV.) 
We fed the centroids of the matched cluster pairs into the IRAF geomap routine, and in 
this way computed second-order polynomial spatial transformation functions. We applied 
the transformations with the IRAF geotran package using bilinear interpolation to rebin 
the images. Due to NGC 6240's drastically different morphology between B and K' bands, 
finding enough clusters in common at these two wavelengths was challenging. Therefore, 
instead of aligning the B image directly to the K' image, the B image was first aligned to the 
I image, and then to the K' image using the transformation computed for the I-K' alignment. 
This two step process also ensures that B — I colors are as accurate as possible, with the 
fewest errors introduced from poor image alignment. Typical rms fitting errors computed by 
geomap were < 0.5 mas when aligning Keck images, < 2.5 mas when aligning HST images, 
and < 5 mas when aligning HST images to Keck images. 

3.4. Aperture Photometry & Flux Calibrations 

Because the star clusters, although unresolved, cannot be conveniently modeled as many 
point sources on a smooth sky background, and because the Strehl ratio for the AO images 
is uncertain, we performed aperture photometry rather than using the standard PSF-fitting 
photometry packages like DAOPHOT. We explored two aperture photometry techniques. 
First, we simply performed aperture photometry on all images using an aperture radius of 
0.07" and a 0.02"- wide sky annulus centered at a radius of 0.1" from the cluster centroid. This 
annulus is near enough to the cluster that we can accurately measure the local background 
in a complicated environment, while the aperture is large enough to lessen the effects of 
poor image alignment. Second, we tried to minimize the color errors by creating matched- 
resolution image pairs before performing aperture photometry. We found that this second 
technique produced poorer results because we were forced to degrade the image quality of 
the NIRC2 images, making many of the clusters disappear into the noise. Therefore, the 
results presented in this paper were calculated by performing simple aperture photometry 
on all the clusters, using the same aperture and annulus size for each filter. Note that the 
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WFPC2 and NIRC2 images have very similar resolutions, with the best resolution (from the 
H band AO image) being a factor of 2 better than the worst resolution (from the WFPC2 
I band image). However the NICMOS images have poorer resolutions than the WFPC2 
and Keck AO images, so that performing aperture photometry on these images with the 
same size aperture and annulus introduces color errors in the case where different colored 
background flux from the galaxy falls within the aperture. Thus caution must be used when 
interpreting the NICMOS data, especially in regions near the nuclei where the background 
flux from the galaxy is bright. 

In Figures [5] - [7] we present the results of our aperture photometry measurements in the 
form of SEDs for each cluster, and in the following paragraphs we outline the procedure we 
used to create these SEDs for two key cases. Other cases can be deconstructed into variations 
on these two cases. 



Case I: Creating an SED for a Cluster Detected Only By HST 

In order to populate the SED of a cluster that is detected in only the B, I, F110W, 
F160W, and F222M images^, first we perform aperture photometry on the cluster in each 
image (using aperture parameters just described). We calculate appropriate aperture cor- 
rections by repeating this procedure on simulated point sources extracted from the Tinytim 
(ver 6.3) software package. Finally, we convert the (aperture corrected) counts to calibrated 
flux (in Janskys) using the PHOTFLAM and PHOTPLAM keywords in the fits header of 
the WFPC2 images, and the instrumental calibration values determined by the NICMOS 
IDT (see §3. lj) . On an SED graph we represent these measured fluxes in the form of a single 
point at the position of the central wavelength for each filter. 

To calculate the photometry errors shown in Figures [5] - [7J for each cluster we simulated 
100 additional clusters using Tinytim PSFs. The simulated clusters were made to have 
the same SED as that which was measured for the real cluster, and they were placed in 
approximately similar areas of the image as the real cluster. (We defined 3 regions of the 
images - far, close, and intermediate distances from the nuclei.) Next, we performed the 
same aperture photometry procedure on the simulated clusters as we did on the real clusters. 
The resulting distribution of 100 measurements is sensitive to the exact placement of the 
simulated clusters and is not gaussian. Therefore we discard the 10 smallest and largest 



2 A cluster would be detected in only the B, I, F110W, F160W, and F222M images, while not being 
detected in the H and K' images, if it is located outside of the NIRC2 narrow camera field of view. Cluster 
#4 is an example of such a cluster for which the described procedure was performed. 
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measured values, and quote errors as 80% confidence limits from the pared distribution. The 
error bars on the single data point at which all models are normalized represent something 
slightly different from the error bars on the other points in the SED. The former are 80% 
confidence limits on the flux at that frequency, while the latter are 80% confidence limits 
on the color. Thus when analyzing these normalized SEDs, which are plotted in log-linear 
space, realize that the SEDs may shift vertically by an amount governed by the error bar on 
the point at which the SEDs are normalized, while the slopes of the SEDs may change only 
by the size of the error bars on the other points. Since the slope of an SED is governed by 
the stellar population of that cluster and any intervening dust, the color error bars indicate 
how confidently we know the cluster's stellar population. Likewise, since the total mass of 
the cluster corresponds to an overall multiplicative factor, i.e. a vertical shift in log-linear 
space, the flux error bars indicate how confidently we know the cluster's mass, assuming no 
errors on the stellar population and extinction values. 

In the error analysis procedure described above, if the 80% confidence limits indicate 
that the cluster may have negative flux, as is sometimes the case when performing aperture 
photometry on barely detected clusters, then we replace the lower error bar with a downward 
arrow, signifying an upper limit. (We still plot the upper error bar.) For clusters that we 
deem undetected we calculate upper limits using a different method. We create a simulated 
cluster using a Tinytim PSF, and place it exactly at the location of the known cluster. We 
decrease the flux of the simulated cluster until it is undetectable (by eye) and call the final 
flux an upper limit. 



Case II: Creating an SED for a Cluster Detected Only with Keck AO 

Because absolute flux calibration of NGS AO images is very difficult when the off-axis 
Strehl is unknown, the procedure we have adopted for populating a cluster's SED with 
photometric information taken by the Keck AO system is quite different from the procedure 
outlined for HST in Case I. As an example, we outline the procedure used to populate the 
SED of a cluster detected only in the H and K' images. Cluster 18 is an example of such a 
cluster for which the following procedure was employed. 

We begin as before by first performing aperture photometry on the cluster in the K' 
image. Next, because we have no way of performing an accurate aperture correction or 
converting counts into Janskys for this image taken using AO, we make use of the NICMOS 
F222M image. Though the particular cluster in question is not visible in the lower resolution 
NICMOS image, one cluster, namely cluster #2, is visible in both images, and is relatively 
far from the bright nuclei. We calculate the total flux of cluster 7^2 in the F222M image 
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using the standard procedure for NICMOS data, namely performing aperture photometry 
and applying an aperture correction calculated from a Tinytim PSF, and converting aperture 
corrected counts into calibrated Janskys by using the appropriate instrumental calibration 
value determined by the NICMOS IDT. Next we perform aperture photometry on cluster 
#2 in the K' image. The ratio of the total flux of cluster #2 to the aperture photometry 
result can be used to convert from counts in the K' image (calculated using a particular 
aperture and annulus) to total, calibrated flux for any cluster. We use this ratio to convert 
the counts in the particular aperture and annulus for cluster #18 in the K' image to total, 
calibrated flux. On an SED graph, this flux is represented in the form of a single point at 
the position of the central wavelength for the K' filter. 

Because the K' and F222M filters are not identical (as can be seen in Fig. [TJ, the 
bootstrapping procedure described above produces accurate results only when the slopes of 
the SED in the NIR are flat, or the same, for clusters 2 and 18. Indeed, the NIR SED is 
very flat for the case of young, dusty clusters . To estimate the e rrors introduced by this 



bootstrapping procedure we convolved various iBruzual fc Charlotl (120031 ) stellar population 
models with the transmission functions for the F222M and K' filters. Assuming a Salpeter 
IMF, solar metallicity, and Padova evolutionary tracks, we determine that observations taken 
with the two different filters should differ by less than 0.1 magnitudes for a cluster between 
the ages of and 10 Gyr. If an intervening dust screen, corresponding to extinctions ranging 
from A v = 4 to 12, reddens the cluster and flattens the NIR SED, then observations taken 
with the two different filters should differ by less than K' AB = 0.05 for a cluster that is less 
than 1 Gyr old. A completely analogous procedure, which makes use of cluster #29, is used 
to bypass the aperture correction and flux calibration of the H band AO images; the errors 
introduced by calibrating H band data with F160W images are comparable to the 0.1 mag 
errors quoted above. 

The photometry error calculations and the upper limit calculations in the case where 
the cluster is only detected in the AO images are completely analogous to the procedures 
described in Case I. However, instead of using Tinytim PSFs to simulate clusters, we used 
a theoretical AO PSF. We model the AO PSFs as the sum of two gaussians. One gaussian 
has a FWHM equal to 0.98X/D where A is the observing wavelength and D is the 10m 
diameter of the Keck telescope. The other gaussian has FWHM equal to 0.98A/r , where 
ro is the Fried parameter, which we assume to be 30 cm. We let the on-axis Strehl ratio 
govern the ratio of the total amounts of energy in each gaussian. We recognize that the 
Strehl used for the simulated clusters is certainly over estimated. One of the consequences of 
these assumptions is evident in Figure [TJ where upper limits for cluster #31 are clearly far 
too low. (Larger Strehls make fainter clusters easier to detect.) Therefore when interpreting 
infrared upper limits shown in Figures [5] - Ui the upper limits computed from NICMOS data 
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- for which the PSF is well modeled by Tinytim - should be trusted more than the upper 
limits computed from the AO data. 



4. RESULTS 

In this complicated galaxy merger the values we have calculated, namely a cluster's 
magnitude in at most 5 distinct wavebands^ are outnumbered by the unknown variables, 
which include, but are not limited to, the following list: 1) the metallicity of gas that 
formed each cluster; 2) the duration and temporal profile of the starburst for each cluster; 
3) the initial mass function (IMF) of stars in each cluster; 4) the stellar evolutionary track 
followed by every star in each cluster; 5) the extinction, scattering, and absorption properties 
of the dust in and in front of each cluster; and 6) the mass and age of each cluster. In 
principle we would like to find the stellar population that most closely matches the observed 
photometry for each cluster, but because it is impossible to fully constrain this problem, 
rather than perform minimization techniques on a many-dimensional parameter space we 
adopt a different approach which requires making the following simplifying assumptions. 
We assume that all clusters: 1) have solar metallicity; 2) were formed in an instantaneous 
burst; 3) are composed of stars with stellar masses governed by a Salpeter IMF; 4) are 
composed of stars that evolve on Padova (1994) evolutionary tracks; and 5) have evacuated 
all intracluster dust before an asre of ^3 and are obscured behind a simple dust screen 

that behaves according to Ry = 3.1. 

The fact that all 32 of these clusters are concentrated within < 5 kpc of the double nuclei 
(in projection) suggests that their formation was triggered by the merger and that they were 
formed roughly contemporaneously. To simplify our analysis we initially assume that all 
clusters were formed simultaneously, although this is obviously false. Later we relax this 
assumption to allow for clusters of three possible young (or intermediate) ages. It is possible 
that some subset of the clusters are old GCs from the progenitor galaxies that happen to be 
passing through the center of the merger, but this is surely not the case for the majority of 
the clusters - especially those closest to the nuclei. Assuming the clusters are not old GCs, 
it is reasonable to choose a metallicity close to solar. If the triggering event was short, then 
the assumption that the clusters' stellar populations are well modeled by an instantaneous 



3 Recall that the F160W and H wavebands are similar, as are the F222M and K' wavebands. 

4 Note that even if some intracluster dust remains after 3 Myr, our assumption is still good assuming that 
the effect of the remaining intracluster dust is far less than that of the dust along the line of sight through 
the merging galaxy. 
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burst is roughly accurate. However, those clusters that are located very near to the double 
nuclei may be experiencing a continual accretion of gas, and therefore it is possible that 
they would be better modeled by a stellar population with a constant SFR. We leave the 
clusters' ages and the amount of intervening dust along the line of sight to each cluster as 
free variables. 

After making the aforementioned simplifications, our approach entails visually compar- 
ing the SEDs that resulted from performing aperture photo metry at B through K' b ands, 



to model SEDs created by the isochrone synthesis code of iBruzual fc Charlotl (120031 1 . In 
Figures [5] - [7| we present the results of this approach for the case where we assume that all 
clusters are 14.5 Myr old. (We justify this choice of age shortly.) We organize the results 
into 3 groups; Figure [5] shows the 12 clusters closest to the northern nucleus, Figure M shows 
the 12 next closest clusters, that lie at an intermediate distance from the norther nucleus, 
and Figure [7] shows the remaining 8 clusters that are farthest from the northern nucleus. 
All model SEDs shown in these figures embody the former 6 simplifying assumptions, and 
for the purpose of visual comparison, each SED has been scaled so that it passes through 
either the K' or the I data pointsjf] The different colors of t he SEP curves correspond to 



varying column densities of interstellar dust applied usi ng the iDrain 3 J2003h d ust pr escrip 



tion with Ry = 3.1. (Grain si ze distribution com es from lWeingartner fc Draind (120011 ); PAH 



optical properties come from lLi fc Draind (120011 ).) For the 12 clusters nearest the northern 



nucleus, column densities corresponding to between 3.5 and 11 magnitudes of visual extinc- 
tion were required to match the observations with the models. Similarly, when comparing 
the intermediate-distance clusters to the model SEDs, column densities corresponding to Ay 
between 4.5 and 14 mag were applied to the models. But for the clusters farthest from the 
northern nucleus, much smaller column densities, corresponding to Ay = 0.5 — 5 mag, were 
needed to produce SEDs spanning the range of observations. Immediately we see that if our 
simplifying assumptions are correct, and if all clusters are assumed for the time being to be 
the same age of 14.5 Myr, then on average the outermost clusters are obscured by a thinner 
dust screen than that which obscures the inner clusters. 

We repeat this procedure for two additional cases, where all clusters are 2.88 Myr old, 
and where all clusters are 180 Myr old. These ages were chosen for the following two reasons. 
1) The ages span a range that is consistent with our previous assumptions. For example, 
clusters younger than 3 Myr could have considerable intracluster dust, and clusters older than 
several hundred million years would have been formed before the onset of the merger, thus 
it is unlikely that they would have formed contemporaneously or be situated coincidentally 



5 All models are for zero redshift. The small rcdshift of NGC 6240 (z = .024) implies an insignificant 
k-correction. 
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near the nuclei 2) The ages bracket the results published in the literature for other star 
clusters in NGC 6240 . Using photometry at the F450W, F547W, and F814W wavebands, 



Pasquali et al. I (120031 ) found that the clusters in NGC 6240's main body and tidal tails have 



probable ages ra nging from 5 t o 13 Myr. By analyzing the CO 2-0 absorption bandhead 



equivalent width, iTecza et al. I (120001 ) show that NGC 6240's K-band light is dominated by 



a ~20 Myr-old, short-duration starburst. 

For each of the three ages analyzed, we choose as the "best-fitting" model the SED with 
the smallest value of Ay that matches the observations for each cluster. For example, the 
observed SED for cluster #13 in Figure [5] is consistent with model SEDs having Ay ranging 
from 5.5 to 16 mag; therefore we assume that cluster #13, if it is 15 Myr old, is extinguished 
by Ay = 5.5 mag. Note that larger extinctions translate to more massive clusters, so our 
method finds the least massive cluster that matches the data for a given age. 

We do not show graphs of the 2.88 and 180 Myr-old model SEDs for every cluster; 
instead we tabulate the results of this exercise in Table [3j The values of Ay listed for 
each cluster and for each assumed cluster age correspond to the "best-fitting" model just 
described. The masses of each cluster are directly proportional to their un-extincted fluxes; 
for two clusters with equal ages and equal apparent magnitudes, the more extinguished 
cluster is more massive. In Figure [8] we depict these same results in the form of histograms. 
The blue, red, and yellow histograms show the "best-fitting" Ay, and corresponding mass 
and (un-extinguished) My, of all 32 clusters for the cases where all clusters are assumed to 
be 2.88, 14.5, and 180 Myr old, respectively. 

While the values listed in Table [3] indicate the stellar population with the "best-fitting" 
Ay assuming 3 different cluster ages, we have not yet specified which age stellar population 
best matches the observations. We address this next step using Figure El in which we use 3 
typical clusters to directly compare the model SEDs with the best-fitting A v for 4 different 
ages. Each of the measured SEDs for clusters #1, 28, and 21 are shown over-plotted with 
model SEDs corresponding to stellar populations that are 2.88, 14.5, 180, and 2000 Myr old. 
Cluster #1 is typical of a cluster for which it is difficult to distinguish between the 14.5 and 
180 Myr-old populations, but it is easy to rule out the youngest and oldest ages with the help 
of spectral information at wavelengths beyond 1/zm. In cases like these, to be consistent with 
our previous method of quoting the least massive clusters that match the data, we choose as 
the most likely cluster age that which yields the smallest cluster mass. For cluster #1, the 



6 Note that in principle the clusters could have been formed during an initial tidal interaction just a few 
hundred million years ago, making them now older than 180 Myr. However, as will be discussed in SJSJ the 
masses required to match the observations to the models for clusters older than ~180 Myr are too large to 
be credible. 
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14.5 and 180 Myr-old populations correspond to ~ 1 x 10 6 M Q and 8 x 1O 6 M , respectively, 
and so we designate this cluster as having aged 14.5 Myr. In Table |3] we indicate this final 
designation with bold text. 

The center panel in Figure [9] uses cluster #28 to demonstrate the case where the 2.88 
Myr-old model SED visibly fits the observations better than stellar populations of the other 
three ages. In cases like these, although the cluster's young age requires it to be very massive, 
we nevertheless designate it as having aged 2.88 Myr by using bold text in Table El Finally, 
cluster #21 is representative of the majority of the clusters, for which it is not possible to rule 
out any age stellar population through visual inspection of the various model SEDs. In cases 
like these, we again resort to choosing the cluster age which yields the smallest cluster mass. 
For cluster #21, the 2.88, 14.5, 180, and 2000 Myr-old populations correspond to clusters 
that have approximate masses of 70 x 1O 6 M , 6 x 10 6 M Q , 60 x 1O 6 M , and 150 x 10 6 M Q , 
respectively, and therefore we designate this cluster (and clusters like this one) as having aged 
14.5 Myr. Note that in order for the younger models to be consistent with the observations 
we must assume more intervening dust; the 2.88 Myr-old model shown for cluster #21 
corresponds to Ay = 11 mag, while the 2000 Myr-old model corresponds to Ay = 6 mag. 
This, of course, demonstrates the typical astronomer's dilemma of distinguishing between 
young, dusty populations and older stellar populations. 

Finally, having just described the method we used to estimate the properties of each 
cluster, we reiterate our motivation to not perform a detailed y 2 minimization when matching 
model SEDs to our observations. Given the number of assumptions we are required to make 
in this work, a detailed minimization of the SED fits would provide very little additional 
insight, and could be misleading. In particular, a minimization would only describe which of 
the three chosen ages best fits the models; a fourth (intermediate) age could fit the models 
better than any of the three. Therefore it is more instructive to examine the cases of young 
and old clusters and to discuss how our results would differ if we had chosen different ages. 
We do this in §5.11 and indeed, we discover - not unexpectedly - that a different intermediate 
age may best describe some of the clusters. 

5. DISCUSSION 

5.1. Properties of the Cluster Population 

As shown with bold text in Table El the best-fitting model stellar populations and dust 
screens imply that the clusters are consistent with being young (~ 3 — 15 Myr) and having 
masses that range from approximately 7 x 10 5 M Q to 1 x 10 8 M Q . The summed mass of all 
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32 clusters is ~ 4 x 1O 8 M , with the single most massive cluster making up ~ 30% of the 
total mass. The 8 next most massive clusters with (1 < M < 6) x 1O 7 M contribute ~ 60% 
of the total mass. Cluster #15, which is coincident with the site of the northern black hole, 
is the only cluster which we expect may have a non-stellar SED, or may be poorly modeled 
by an instantaneous burst. This cluster is the second most massive cluster and accounts for 
approximately 10% of the total mass. 

With the exception of the most massive cluster, the clusters analyzed in this paper 
have mas ses similar to, albeit sl ightly larger than, those found in other well-studied merging 
galaxies. IWhitmore &: Zhangl (120021 ) estimate that the young, optically bri ght clusters in 



the A ntennae galaxies have masses ranging from 2 x 1O 4 M to 4 x lO e M . iGilbert et al. 
(120001 ) find that the brightest near-IR star cluster in the Antennae is ~ 4 Myr old and has 
M = 1.6 x 1O 7 M . Similarly, many of the young (< 10 Myr) c lusters in Arp 220 have 
masses > 1O 6 M , or even as much as 1O 7 M ( Wilson et al. 12006 ). (Note that all three of 
these aforementioned works assume a Salpeter IMF.) iMcCrady et al.1 (120031 ) found kinematic 
masses of 3.5 x 1O 5 M and 1.5 x 1O 6 M for two clusters in M82. 

Recall that the masses stated in this paper are underestimated in the sense that we 
adopted stellar population models with the smallest value of Ay possible. In addition, for 
all cases where the observed SEDs fit the different aged models equally well, we designated 
the correct model to be that which yielded the least massive cluster. But by comparing the 
resulting cluster masses to those in the literature, we judge that these two assumptions were 
probably necessary for the majority of the clusters. Figure [8] shows that if we had assumed 
the clusters were much younger or older than ~ 15 Myr, many of the resulting masses would 
be far larger than those measured in any other merging system. In fact, the clusters' masses 
put stringent constraints on their maximum ages, assuming we have correctly modeled other 
properties of the stellar populations. As clusters age past 180 Myr, they become dimmer 
and dimmer so that the masses required for us to detect them become unreasonably large 
(~ 1O 9 M ). 

The star formation rates measured for NGC 6 240 using a var i ety of methods put ad- 
ditional stringent constraints on the cluster's ages. iBeswick et al. I (120011 ) use the 1.4-GHz 
lumin osity to show tha t _SFR(M > 1M ) = 83.1 M yr _1 in the nuclear region of NGC 
6240. iGao fe Solomon! (120041 ) show that the HCN luminosity, which traces dense molecu- 
lar gas, is a good indicator of the star format i on rat e in spirals and ULIRGs alike. Using 
the relationship descr ibed by iGao fc Solomon! (12004 ) and the value of L H cn measured by 
Solomon et al. I (119921 ). we calculate that the total star formation rate in NGC 6240 is ap- 
proximately 200 M yr _1 . This value is sl ightly larger than tha t derived from NGC 6240's 



integrated IR luminosity, ~ 140 M yr 1 (IHeckman et al. Ill990l ) 
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The total star formation rate due to the 32 clusters analyzed in this paper is ~ 9 yr 1 



Considering that our measurement exceeds the SFR measured by iBeswick et al. I ( 120011 ) for 
the nuclear region of NGC 6240, and that the unresolved clusters do not include the major- 
ity of the infrared nuclear light (see Fig. [2D, our measured SFR may be implausibly high. 
However, ~ 80% of the total SFR from the clusters comes from just 4 clusters which each 
contribute SFR > 6 M Q yr -1 . These 4 clusters are the only clusters for which we determined 
that the 3 Myr-old model SEDs fit the observations better than models of other agesQ We 
have already conceded that our method of visual inspection does not search the entire pa- 
rameter space for a best fitting model. While a Salpeter IMF or a solar metallicity may be to 
blame, it is unlikely that all but 4 clusters can be reasonably modeled by our assumptions. 
Thus we suggest that these 4 clusters are neither 3, nor 15, nor 180 Myr old, but rather 
some other intermediate age. The K band brightn ess implied by the instantaneous burst 



stellar population models (IBruzual &: Charlotl 120031 ) monotonically decreases for population 



ages greater than ~ 10 Myr, and increases for population ages between ~ 3 and 10 Myr. 
Thus, these 4 clusters may be better described by a stellar population with an age of ~ 10 
Myr, implying that they are less massive and contribute less to the total SFR. 

The total SFR measured for the remaining 28 clusters, for which the observations are 
consistent with 15 Myr-old stellar populations, is ~ 15M Q yr -1 . This agrees well with the 
literature, and implies that the newly detected clusters account for ~ 20% of the nuclear 
SFR, or ~ 10% of the total SFR in the ULIRG. Finally, if we assumed that the same 28 
clusters were all 3 Myr old and summed their appropriate star formation rates in Table [31 we 
would conclude that they contribute a whopping ~ 879M yr -1 . Thus given our assumptions 
(outlined in detail in paragraph 1 in §3]) the total star formation rate for NGC 6240 puts 
rather stringent constraints on the clusters' ages. 

To summarize, the masses of star clusters observed in other mergers like Arp 220 and 
the Antennae strongly suggest that the clusters in NGC 6240 are younger than 180 Myr, 
while the SFRs measured for NGC 6240 using Li. 4 gh z , Lhcn, and L 8 _ioooAtm require that the 
majority of the clusters be older than 3 Myr. 28 of the 32 detected clusters are consistent with 
being 15 Myr old, while the large SFRs of 4 clusters (clusters #6, 28, 29, and 32@) suggest 
that these 4 may be more consistent with another intermediate age. The large mass of cluster 
#15 and the fact that cluster #15 is coincident with the northern black hole suggests that 
the measured SED may have a non-stellar component. NGC 6240's cluster population - 



7 A11 four of the of the clusters that have SFR > 6 M 



also have M > 10 7 M 



©• 



3 Note that because clusters #29 and 32 lie outside of the "nuclear" region for which lBeswick et al. ( 2001) 



measu red L 14 GHz , their high SFRs do not directly contradict 83.1 M Q yr 1 limit set by 
(l200ll) . 



Beswick et al. 
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excluding the 5 aforementioned clusters - is very similar to other cluster populations in the 
literature. The cluster masses range from approximately 7 x 1O 5 M to 4 x 10 7 M Q , and the 



total contribution to the SFR from these 27 clusters is 
SFR in the nuclear region. 



10M Q yr -1 , or ~ 10% of the total 



The cluster populat ion discussed in this p aper is considerably more massive than that 
which was discussed in iPasquali et al. I (120031 ). The latter consists of 41 clusters located 



i n NG C 6240's "main body" and 13 clusters located in the galactic tails. IPasquali et al. 
J2003TI estimate that these clusters have probable masses of (1 — 2) x 10 5 M e . Because 



Pasquali et al. ( 2003 ) do not give the locations of the clusters they analyz^l it is difficult 
to know whether this is a discrepancy or whether we ar e simply noticing an in crease in the 
cluster mass toward the double nuc lei, as was noted by Pasquali et al. I (j2003r ). While it is 
obvious that IPasquali et al. I (120031 ) did not analyze the majority of the clusters discussed 
in this paper, for which we measure only upper limits on the B and I photometry, it is 
intriguiging that clusters #2 and 4, which are located in what appears to be the "main 
body" of NGC 6240, and for which B and I measurements were m ade, are rather massive 
(> 10 6 M ) compared to estimates made b ylPasauali et al. I (|2003l ) . It is possible that this 
apparent discrepancy is due to the fact that Pasquali et al. ( 20031 ) used Starburst99 stellar 
population models with a different reddening law and sub-solar metallicity. This apparent 



discre pancy cannot, however, be attributed to a difference in IMF, since IPasquali et al. 



(120031 ) used the same IMF as we did in this paper, namely a Salpeter IMF with lower and 
upper mass cutoffs of 0.1 and 100 M , respecitvely. 

Finally, we note that wh il e usin g a Salpeter IMF is helpful when comparing our results 
to the literature, the iKroupal (120011 ) IMF is more realistic and results in lower cluster mass 
estimates. Since the infrared luminosity of clusters is dominated by light from high mass 
stars that quickly turn off the main sequence, and since the Salpeter and lKroupal (120011 ) IMFs 
differ only for masses below 0.5 M , the cluster masses one derives when employing these 
IMFs differ by a simple scaling factor for clusters with ages less than the main sequence life- 
time of a 0.5 M Q star. The exact mass scaling can be easily derived by piecewise-integrating 
the two IMFs from 0.1 to 100 M Q . The result is that two young clusters with equal infrared 
luminosities but different IMFs have masses that differ by a factor of ~ 1.3, with the larger 
mass corresponding to the cluster forme d with a Salpe ter IMF and the smaller mass cor- 
responding to the cluster formed with a IKroupal (120011 ) IMF. Thus if we had used a more 
realistic IMF throughout our analysis in this paper, we would have estimated cluster masses 



Pasquali et all |2003h did not publish the locations of the clusters they analyze, nor did they provide 



this information on the website they mention. We were not able to obtain the cluster positions from private 
communication with the authors. 
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to be smaller by a factor of about 0.77. Correspondingly, our SFR estimates would decrease 
by this same factor. 



5.2. Dust and Extinction in NGC 6240 



NGC 6240' s cluster popu lation, if in fact obscured by a simple dust screen with Ry = 3.1 
as described by iDraind (l2003h . is extinguished by column densities of dust corresponding to 
Ay between 0.5 and 14 mag There is no significant difference in the values of Ay measured 
for the 24 clusters closest to the northern nucleus, but the 8 clusters farthest from the nucleus 
exhibit significantly less extinction. This is to be expected for an active merger, in which 
gas is tunneled onto the nuclei, in this case powering the AGN that iKomossa et al. I (120031 ) 
detected. Excluding clusters #6, 28, 29, 32, and 15 as previously discussed, on average the 
6 clusters farthest from the northern nucleus are extinguished by Ay = 2.3 mag, and the 
remaining 21 clusters closest to the nucleus are extinguished by A v = 7 mag. 



Gerssen et al. I ( 120041 ) and iTecza et al. I (120001 ) both estimate extinctions toward the 



northern and southern nuclei assuming foreground dust screens. For the northern nucleus 
we can directly comp are our results to t heirs using clu sters #12 and 14 as two sight lines 
through this region. (IGerssen et al. I and ITecza et al. I both class ify the brightest northern 



Gerssen et al. 



and 



Tecza et al. estimate that A 



v 



2.35 



knot as the northern nucleus.) 
mag and Ay = 1.6 mag, respectively. Both of these values are considerably less than the 
measurements made in this paper: Ay = 6.5 ma g and Ay = 3.5 mag toward clusters #12 
and 14, respectively. In addition, ITecza et al. I measure a peak extinction of Ay = 7.2 
mag, while we find that man y clusters discu ssed in this pape r have Ay > 7 mag. For 



these reasons we suggest that IGerssen et al. I and ITecza et al. I have both underestimated 
the amount of extinction in this ULIRG. Our values of extinc tion agree bet t er wit h the 



values in the literature, measured using a variety of methods. iVignati et al. I (119991 ) used 
X-ray observations and found Ay ~ 10 mag, iGenzel et al. I (119981 ) used line ratios and found 
Ay > 5 mag, and iRieke et al. I (119851 ) used continuum K and L spectrophotometry and 
inferred Ay ~ 15 mag. The cause of the previous underestimates may be due in part to 
the different spatial scales being probed, and in part to the fact that the measurements can 
only probe the optically thin layer surrounding the nuclei. The observations presented in 
this paper are the first infrared observations to resolve the northern nucleus, and thus it 



is sensible that we would infer higher values of extinction than were found by ITecza et al. 



10 We have not accounted for Galactic extinction when estimating these values. The Galactic extinction 
toward NGC 6240 is Ay = 0.25 mag and is negligible given the uncertainty on our measurements. 
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feoooh . 

Finally, note that while iTecza et al. I (120001 ) may have underestimated the overall extinc- 
tion in this ULIRG, the morphology of their extinction map is consistent with our findings. 
They find that the extinction is greatest between the nuclei and least toward the north- 
ern nucleus. We detect a cluster coincident with the northern nucleus and we find a large 
concentration of clusters just northwest of the northern nucleus, yet we cannot resolve the 
southern nucleus, nor do we detect many clusters directly between the nuclei. Thus the 
locations of the clusters in our sample strongly suggest that the region between the nuclei is 
the most obscured, and that the southern nucleus is more highly obscured than the northern 
nucleus. This evidence suggests that deeper observations would result in additional cluster 
detections, particularly between the two nuclei. 



5.3. Summary 

The circumnuclear locations of the clusters detected in these new, high- resolution obser- 
vations strongly suggest that their formation was triggered by the merger event that created 
the ULIRG NGC 6240. This motivates our decision to model the clusters as solar metallicity, 
instantaneous burst stellar populations. If we further assume that the stellar populations 
were created with a Salpeter IMF, and that the stars in the clusters follow Padova evolution- 
ary tracks, then we find that 27 of the 32 newly discovered clusters are consistent with being 
~ 15 Myr old, and 4 of the 32 clusters are probably between 3 and 15 Myr old. Cluster ages 
much younger or much older than 15 Myr can be ruled out by the large star formation rates 
and masses inferred, respectively. Assuming the clusters are 15 Myr old, then the majority 
of the clusters have masses ranging from 7 x 1O 5 M to 4 x 1O 7 M , and the total contribution 
to the SFR from these 27 clusters is ~ 10M Q yr _1 , or roughly 10% of the total SFR in the 
nuclear region. 

The extinctions calculated toward these clusters range from Ay = 0.5 - 14 mag, with the 
least extinguished clusters being located farthest from the double nuclei. While the range 
of extinctions quoted in the literature for NGC 6240 varies greatly, the values measured 
in this paper agree well with some of the larger published values. The locations of th e 



clusters detected here, along with the extinction map published by ITecza et al. I (120001 ). 
suggest that there may be as- yet-unobserved clusters located between the nuclei that are still 
obscured. The complicated morphology of this ULIRG makes it very difficult to determine 
our completeness limits empirically. 
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Table 1. Keck AO Observations of NGC 6240 using NIRC-2 Narrow Camera 



Date 


Filter 


Integration Time (s) 


Guide Star's 








On-Axis Strehl 


24 Apr 2005 


H 


800 


0.21 


17 Aug 2003 


K' 


630 


0.32 



Table 2. Archival HST Observations of NGC 6240 



Observing Program 


Instrument 


Camera 


Filter 


Exposure Time (s) 


GO 6430 (PI: van der Marel) 


WFPC2 


PC 


F814W (B) 


3 x 400 




WFPC2 


PC 


F450W (I) 


3 x 700 


GTO 7219 (PI: Scoville) 


NICMOS 


NIC2 


F110W 


4 x 40 




NICMOS 


NIC2 


F160W 


4 x 48 




NICMOS 


NIC2 


F222M 


4 x 56 



Table 3. Identification of Infrared Clusters in NGC 6240 



Cluster 
ID # 


A X 

( E) 


A Y 

( N) 


M v (mag) 
2.9, 14.5, 180 Myr 


A v (mag) 
2.9, 14.5, 180 Myr 


M/(10 D M©) 
2.9, 14.5, 180 Myr 


SFR (M yr : ) 
2.9, 14. 5, 180 Myr 


1 

X 


1 43 


9 95 

U . UtJ 


-14.9, 


-13.0, 


-13.6 


2 5 5 10 

Z;.(J) VJ •KJ . X.W 


2 8 9 a 7 9 


97 06 04 

U.J 1 j VJ •VJVJ . W.Wt: 


2 


-0 96 

W . U\J 


4 89 


-18.0, 


-15.8, 


-16.4 


4 5 2 2 5 


49 1 1 4 a 1 05 7 

t: <_/ . W . J — L • t: ^ XWW. 1 


1 7 00 79 59 

x i .ww^ vj » i kJ • yj.vJij 




-0 06 

W . WW 


2 23 

Z. . Z. O 


-17.6, 


-14.1, 


-14.8 


12 6 5 7 


S4 2 4 b 23 9 


1 1 81 17 13 

XX.VJXj VJ • X 1 . U.1U 


4 


-7 31 

1 .OX 


1 87 

X . VJ 1 


-16.7, 


-14.5, 


-15.1 


6 3 5 4 


14 8 3 4 a 31 9 


5 13 24 18 

a.lUj W • Z^ t: . W.XVJ 


W 


01 

W . W X 


1 34 

X . U1 


-19.1, 


-15.8, 


-16.1 


18 13 5 110 


137 5 11 1 b 77 3 


47 73 77 43 


6 


0.41 


1.12 


-17.7, 


, -14.4, 


-15.1 


9 5 6 

%J • V . u • w « w.w 


35 8 3 31 4 

*JVJ • VJ ' - O • W . 'J -L • i 


12 42 21 17 


7 


-0 80 

w . ow 


63 

W . WO 


-17.2, 


-13.8, 


-14.3 


10 5 6 5 

lU.Uj VJ »\J . xj .W 


22 9 1 8 b 15 2 

Z. Z. . .J . X . O . X U . Z. 


7 95 13 08 

I .tyvJ^ W • X O . W.WVj 


a 


-0 15 

W . X W 


44 

w . n 


-15.7, 


-12.4, 


-13.1 


7 3 5 4 

i.W. ti.W 


5 6 5 b 4 9 


1 93 03 03 

X.<_/0^ W • W »J • w.wo 


Q 


-0 46 

W . t: W 


0.42 


-18.0, 


-14.6, 


-15.4 


10 5 6 7 5 

XW.W. VJ • VJ • t • ') 


48 8 3 9 b 42 8 


1 6 93 27 24 

1 V . -J . VJ • Ad 1 . W . Zjt: 


10 

X W 


-0 67 

W . W 1 


41 

W . t: X 


-16.7, 


-13.3, 


-14.0 


9 5 5 5 5 

<_/.W^ %J»\Jy O.O 


13 8 1 1 b 10 8 

XO.VJ^ XaX . XW.VJ 


4 80 08 06 


11 


-0 82 


32 

W . (J Zj 


-18.3, 


-15.0, 


-15.4 


15 11 9 


63 5 5 4 b 39 9 

V - J . / . iJ • j: . - J : / . • / 


22 05 37 22 


1 2 


26 

W . Zj W 


17 

W. X 1 


-20.4, 


-17.1, 


-17.8 


10 6 5 6 5 

XW.Wj VJ • <J * W . "vJ 


433 5 39 b 359 6 

t:*JO . Wj Ui/<U . OUl/.U 


1 50 53 2 70 1 99 


1 3 

X *J 


-0 25 

W . ZjO 


1 5 

W. XW 


-17.6, 


-14.3, 


-14.9 


10 5 5 6 

XW.Wj o • o . w.w 


34 6 2 8 b 27 2 


1 2 02 19 15 

IZj.UZj, W • X iJ . W.XW 


14 


28 

W . Zj VJ 


11 

w . x x 


-20.1, 


-16.7, 


-17.4 


8 3 5 4 

VJ.W^ *J • Kj . T: . W 


325 2 26 2 b 255 2 


1 1 2 93 181 141 


J. U 


00 

W . W W 


00 

W . WW 


-21.0, 


-17.5, 


-18.2 


14 5 9 9 5 


782 1 56 3 b 549 1 


271 57 3 90 3 04 


16 


-0.41 


-0.06 


-19.1, 


-15.5, 


-16.2 


14.0, 8.0, 8.5 


129.7, 8.8 b , 86.2 


45.05, 0.61, 0.48 


17 


-0.27 


-0.13 


-18.6, 


-15.1, 


-15.9 


12.0, 6.5, 7.5 


86.2, 6.2\ 64.0 


29.93, 0.43, 0.35 


18 


-0.53 


-0.37 


-18.9, 


-15.2, 


-16.0 


15.0, 8.5, 10.0 


105.9, 6.8 b , 74.3 


36.77, 0.47, 0.41 


19 


0.06 


-0.94 


-16.8, 


-13.4, 


-14.1 


11.0, 6.0, 6.5 


16.4, 1.3\ 12.2 


5.71, 0.09, 0.07 


20 


0.19 


-0.97 


-17.3, 


-13.9, 


-14.6 


12.5, 7.5, 8.0 


26.1, 2.0 b , 19.4 


9.06, 0.14, 0.11 


21 


-0.45 


-1.15 


-18.4, 


-15.1, 


-15.8 


11.0, 6.5, 7.5 


71.6, 5.8\ 59.4 


24.85, 0.40, 0.33 


22 


-0.62 


-1.19 


-16.2, 


-13.0, 


-13.7 


8.5, 5.0, 6.0 


9.3, 0.8\ 8.7 


3.24, 0.06, 0.05 



Table 3 — Continued 



Cluster 


A X 


A Y 


M v (mag 


) 


A v (mag) 


M/(1O 6 M ) 


SFR (M yr- 1 ) 


ID # 


("E) 


("N) 


2.9, 14.5, 180 Myr 


2.9, 14.5, 180 Myr 


2.9, 14.5, 180 Myr 


2.9, 14.5, 180 Myr 


ZO 


-l.U ( 


-i.zy 


-17.6, -14.3, 


-15.0 


1 1 n t n i n 
11. U, /.U, l.U 


Q/l 1 O Ob 07 O 
o4. ( , Z.y , Z I .Z 


in n on n i r 
lz.U4, U.ZU, U.lO 


z4 


l.»4 


i fin 

-l.OU 


-14.4, -12.5, 


-12.8 


/in on on 
4.U, Z.U, z.U 


ion K^> Q fi 
l.o, U.D , o.o 


n si n r\A n no 
U.Ol, U.U4, U.Uz 


25 


0.13 


-2.21 


-16.7, -13.5, 


-14.2 


8.0, 5.0, 5.5 


14.5, 1.4 b , 13.4 


5.03, 0.10, 0.07 


26 


-0.22 


-2.46 


-15.9, -12.7, 


-13.2 


8.5, 5.5, 4.5 


7.2, 0.7\ 5.7 


2.51, 0.05, 0.03 


27 


0.89 


-2.62 


-18.0, -14.6, 


-15.3 


19.0, 14.0, 15.0 


50.0, 3.8\ 39.2 


17.35, 0.26, 0.22 


28 


-0.19 


-3.10 


-17.9, -14.6, 


-15.3 


8.0, 4.5, 5.5 


42.3, 3.8, 39.2 


14.68, 0.26, 0.22 


29 


6.58 


-4.13 


-17.0, -14.5, 


-15.1 


6.5, 3.5, 4.0 


19.3, 3.4, 31.9 


6.71, 0.24, 0.18 


30 


0.95 


-4.25 


-16.0, -13.2, 


-13.8 


7.5, 4.0, 4.5 


7.3, 1.0\ 9.2 


2.53, 0.07, 0.05 


31 


-1.02 


-5.30 


-14.4, -12.8, 


-13.4 


3.5, 2.0, 2.5 


1.8, 0.7 a , 6.6 


0.62, 0.05, 0.04 


32 


-3.54 


-7.71 


-19.0, -16.2, 


-16.9 


7.5, 4.0, 4.5 


124.7, 17.0, 157.6 


43.31, 1.18, 0.87 



Note. - Column 1 specifies each cluster using the identification numbers as shown in Figure HI Smaller 
numbers correspond to more northern clusters. Columns 2 and 3 give the clusters' positions relative to cluster 
#15, which is coincident with the site of the northern black hole. Column 4 shows the unextincted V band 
absolute magnitude each cluster would have if it were 2.9, 14.5, and 180 Myr old, and if it were extincted by 
an amount specified in column 5. (See text for detailed description of the assumptions made for the stellar 
populations.) Column 5 specifies the minimum amount of extinction needed for the measured SEDs to match 
the model SEDs for the 3 cluster ages. Column 6 gives the mass of each cluster for the 3 cluster ages assuming 
the values of Ay listed in column 5. Errors on the mass and extinction can be deduced from the error bars 



shown in Figures [5E Column 7 shows the star formation rate implied by the masses and ages listed in column 
6; we simply divide the mass by the age in each case. Bold text indicates which of the three stellar population 
ages best agrees with our observations, as determined through visual inspection. 

a Using visual inspection of model SEDs, it is difficult to distinguish between the 14.5 and 180 Myr-old stellar 
populations, while we can rule out the 3 Myr population. In these cases, to be consistent with our method 
of quoting the minimum amount of extinction needed for the measured SEDs to match the model SEDs, we 
choose as the most likely cluster age that which yields the smallest cluster masses. The 14.5 Myr-old stellar 
population implies a less massive cluster. 

b Using visual inspection of model SEDs, there is insufficient data to distinguish between stellar populations 
with ages ranging from 2.9 to 180 Myr. In these cases we choose as the most likely cluster age that which yields 
the smallest cluster masses. The 14.5 Myr-old stellar population corresponds to the least massive cluster. 
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Wavelength [jj,m] 

Fig. 1. — Filter transmission functions for the NICMOS (solid line) and Mauna Kea (dotted 
line) filter sets. The curves centered around 1.6 /xm are the F160W and H filters and the 
curves centered around 2.2 /im are the F222M and K' filters. 
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Fig. 2. — Left and middle panels show logarithmically scaled WFPC2 images at B and I 
bands, respectively. Known cluster locations are marked with circles. '+' marks the location 
of cluster #15 which is coincident with the northern black hole. Images are 2.5" x 5" (1.2 x 
2.4 kpc); North is up and East is to the left. Right panel shows the same cluster locations 
superposed on a logarithmically scaled, unsharp masked K' image taken with NIRC2 using 
Keck NGS AO. The northernmost and southernmost clusters encircled are clusters #5 and 
#28, respectively. 



31 




Relotive Arcsec 



Fig. 3. — Relative fields of view for the complete set of 7 images analyzed in this paper. The 
blue square shows the WFPC2 Planetary Camera FOV (F450W and F814W filters). The 
purple square shows the NICMOS NIC2 FOV (F110W, F160W, and F222M filters). The red 
and green squares show the NIRC2 narrow camera FOV (K' and H filters, respectively). Red 
dots show the locations of all 32 clusters analyzed in this paper, for which we have spectral 
information in at least one of the NIR wa vebands. The oran ge cross marks the location of 
the northern black hole, as determined by Max et al. I (120061 ). The overlaid grayscale image 
shows NGC 6240 taken with the WFPC2 F450W filter. 
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o 

ID 




Relotive Arcsec 

Fig. 4. — Position and identification number of all 32 clusters, shown relative to the infrared 
fields of view. Scale in arcseconds is ~ 1/4 that in Fig. [3J Identification numbers increase 
toward the south. The orange cross indicates the position of the northern black hole which 
is coincident with cluster #15. 
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Fig. 5. — Open diamonds, open squares, and filled circles mark the photometric measure- 
ments made for 12 clusters using WFPC2, NICMOS, and NIRC2 images, respectively. The 
12 clusters depicted here are those situated nearest to the northern nucleus, in projection; 
cluster labels are given in the lower right corner of each panel. Error bars show 80% con- 
fidence limits. Upper limits are depicted with a downward arrow for those filters where a 
cluster is not detected. Read §3.41 to understand why some clusters show both upper error 
bars and upper limits for a single filter. Over-plotte d on the measured d a ta po ints are model 
SEDs created by the population synthesis code of iBruzual fc Charlotl (120031 ). The model 
SEDs show 14.5 Myr-old instantaneous burst, solar metallicity, stellar populations, assum- 
ing a Salpeter IMF and Padova (1994) evolutionary tracks. Models hav e been extingu ished 
by varying column densities of dust in a simple dust screen, following the iDraind (120031 ) dust 
prescription. Red, orange, green, and blue curves correspond to Ay = 3.5, 6.5, 8.5, and 11 
mag, respectively All model SEDs, as well as the observed K' photometry, are normalized 
to the K' data point for ease of comparison. 
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Fig. 6. — Similar to Figure [5], but here we show the 12 next closest clusters to the northern 
nucleus. Red, orange, green, and blue curves correspond to 14.5 Myr old stellar populations 
with Ay = 4.5, 6.5, 8.5, and 14 mag, respectively. 
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Fig. 7. — Similar to Figure [5l but here we show the 8 clusters farthest from the northern 
nucleus, with all photometry and model SEDs normalized to the I data point for ease of 
comparison. Red, orange, green, and blue curves correspond to 14.5 Myr old stellar popu- 
lations with A v = 0.5, 2.0, 3.5, and 5 mag, respectively. Read §3.41 (case II) to understand 
why the H and K' upper limits shown for cluster #31 should not be trusted. 
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Fig. 8. — Red histograms show the un-extincted absolute V-band magnitudes, masses, and 
extinctions of clusters when all clusters are forced to be 14.5 Myr-old, and differ only by 
the amount of extinction applied to each. (e.g. Red histograms summarize the results of 
Figures [5H3 for the best-fitting Ay.) Blue and yellow histograms show similar summaries for 
the cases where we force all clusters to be either 2.88 Myr or 180.5 Myr old, respectively. 
Note that these histograms include data for all 32 detected clusters, although we believe 5 
of the clusters may be better described by different, intermediate ages. 
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10.00 




Wavelength (//m) 

Fig. 9. — Open diamonds, open squares, and filled circles mark the photometric measure- 
ments made for clusters #1, 28, and 21 using WFPC2, NICMOS, and NIRC2 images, re- 
spectively. Over-plotted SEDs in red, orange, green, and blue show 2.88, 14.5, 180, and 2000 
Myr-old instantaneous burst stellar populations, respectively. For each model SED plotted, 
we have applied the amount of extinction necessary to make the SEDs most closely match 
the observed photometry. For cluster #1, the red, orange, green, and blue curves correspond 
to A v = 2.5, 0.5, 1, and mag. For cluster #28, the red, orange, green, and blue curves 
correspond to Ay = 8.0, 4.5, 5.5, and 4 mag. For cluster #21, the red, orange, green, and 
blue curves correspond to Ay = 11, 6.5, 7.5, and 6 mag. Read Figure [5] caption for more 
details. 



